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Abstract. We develop the quantum Boltzman equation approach for the Kondo 
breakdown quantum critical point, involved with two bands for conduction electrons 
and localized fermions. Particularly, the role of vertex corrections in transport is 
addressed, crucial for non-Fermi liquid transport of temperature linear dependence. 
Only one band of spinous may be considered for scattering with gauge fluctuations, and 
their associated vertex corrections are introduced in the usual way, where divergence 
of self-energy corrections is cancelled by that of vertex corrections, giving rise to 
the physically meaningful result in the gauge invariant expression for conductivity. 
On the other hand, two bands should be taken into account for scattering with 
hybridization excitations, giving rise to coupled quantum Boltzman equations. We 
find that vertex corrections associated with hybridization fluctuations turn out to be 
irrelevant due to heavy mass of spinous in the so called decoupling limit, consistent 
with the diagrammatic approach showing the non- Fermi liquid transport. 
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1. Introduction 

Non-Fermi liquid transport phenomena near heavy fermion quantum critical points 
(QCPs) are one of the central interests in the business of strongly correlated electrons 
[l]. To reveal scattering mechanism between charge carriers and critical fluctuations is 
the key to understand the temperature linear resistivity [2], one of the hallmarks for 
non- Fermi liquid physics in the quantum critical matters. Basically, the nature of critical 
modes, more concretely, the dynamical exponent z determining their dispersion relation 
uj and the scattering vertex between charge carriers and such critical modes are 
essential ingredients for the transport phenomena. 

The standard model of quantum criticality in a metallic system is a z = 2 critical 
theory, often referred as Hertz- Moriya-Millis (HMM) theory [3]. Fortunately, many 
heavy fermion compounds have been shown not to follow the spin-density-wave (SDW) 
theoretical framework, where the temperature linear resistivity [1], divergent Griineisen 
ratio with an exponent 2/3 [5|, Fermi surface reconstruction at the QCP [HI [7], and the 
presence of localized magnetic moments at the transition towards magnetism [8] seem 
to support a more exotic scenario. An interesting suggestion is that the heavy-fermion 
quantum transition is analogous to an orbital selective Mott transition P, [101 [HI [12], 
where only the f-electrons experience the metal-insulator transition, identified with a 
breakdown of the Kondo effect. This Kondo breakdown scenario differs from the HMM 
theory [3], in respect that the whole heavy Fermi surface is destabilized at the QCP in 
the former case while only hot regions connected by SDW vectors become unstable in 
the latter. 

The nature of the Kondo breakdown QCP turns out to be multi-scale [TTl [12] . 
Dynamics of hybridization fluctuations is described hj z = 3 critical theory due to 
Landau damping of electron-spinon polarization above an intrinsic energy scale E*, 
while by z = 2 dilute Bose gas model below E*. The energy scale E* originates from 
mismatch of the Fermi surfaces of conduction electrons and spinous, shown to vary from 
(9(10°) mK to 0(10^) mK. Based on the z = 3 quantum criticality, both the logarithmic 
divergent specific heat coefficient and power-law diverging thermal expansion coefficient 
was shown to fit successfully, giving rise to the divergent Griineisen ratio with an 
exponent 2/3 [13] . 

Recently, we showed that not only electrical resistivity but also thermal resistivity 
shows quasi-linear temperature dependence around the Kondo breakdown QCP due 
to scattering with z = 3 critical hybridization fluctuations [H], based on the Kubo 
formula where diagrammatic calculations were performed in the 1/A^ expansion with 
the spin degeneracy A^. An important point is that vertex corrections for scattering 
with hybridization fluctuations can be neglected, a unique feature of the two band 
model, resulting from heavy mass of spinous [TTl [T2] . This allows us to replace the 
transport time with the scattering time for such a process. Then, only the self-energy 
correction will be involved for transport, causing the non- Fermi liquid resistivity. 

In this paper we clarify the issue related with vertex corrections for the transport 
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phenomena at the Kondo breakdown QCP, based on quantum Boltzman equations 
where vertex corrections are introduced naturally. An important feature is emergence of 
coupled quantum Boltzman equations for the distribution function of each band. In the 
heavy fermion phase described by condensation of Kondo bosons, only the lowest heavy 
fermion band may be taken into account, resulting in the Fermi liquid transport owing 
to the absence of scattering with gapless fluctuations. In the fractionalized Fermi liquid 
phase of the Kondo breakdown scenario [10], hybridization fluctuations are gapped, 
leaving the two bands decoupled in the low energy limit and allowing us to consider 
the two quantum Boltzman equations independently. On the other hand, at the Kondo 
breakdown QCP critical hybridization fluctuations force the two quantum Boltzman 
equations coupled, requiring to take both distribution functions int account on equal 
footing. In this situation the quantum Boltzman equation study was not preformed yet 
at least for the heavy fermion QCP as far as we know. 

The present paper is organized as follows. In section 2 we introduce an effective 
field theory for the Kondo breakdown QCP. In section 3 we examine the electrical 
resistivity based on the coupled quantum Boltzman equations, where the quantum 
Boltzman equation study is reviewed for the U(l) gauge theory of one band model 
and its extension to the two band model is derived. In section 4 we summarize our 
results. 

2. An effective field theory for the Kondo breakdown quantum critical 
point 

We start from the Anderson lattice model in the large-?/ limit 

« (ij) 

+ VY.{dlc,^ + H.c.) 

i 

+ dlid^ + ef)d,^ + Jj2Si- S„ (1) 

« (ij) 

where Qo- and di^r are conduction electron with a chemical potential fi and localized 
electron with an energy level ef. The last spin-exchange term is introduced for 
competition with the hybridization term of V. 

Resorting to the U(l) slave-boson representation di„ = h\fi„ with the single 
occupancy constraint h\hi + fl^fia = SN to take strong correlations with S = 1/2, 
one can rewrite Eq. (1) with 

L = Y1 4a(^r - f^)Cia " tY,{ctcj^ + H.C.) 

+ VY.ibiflQa + H.C.) + Y.bldrk 

i i 

+ E flidr + ef)U + ^ E(/ilx.,/,. + H.C.) 
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+ 1 y: Hbik + flu - 1) + ^ E \x^J\^ (2) 
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where bi and fia are holon and spinon, associated with hybridization and spin 
fluctuations, respectively. The spin-exchange term for the localized orbital is 
decomposed via exchange hopping processes of spinous, where Xij is a hopping parameter 
for the decomposition. Aj is a Lagrange multiplier field to impose the constraint, and 
N is the number of fermion flavors with a = 1, A^. 

Performing the saddle-point approximation of hi b, Xij ~^ Xi and iXi A, one 
finds an orbital selective Mott transition as breakdown of the Kondo effect at J ~ Tx, 
where a spin-liquid Mott insulator {{bi) = 0) arises in J > Tk while a heavy-fermion 
Fermi hquid {{k) ^ 0) resuhs in Tf^ > J [TOl HH [12]. Here, Tk = Dexp(^^^) is 
the single impurity Kondo temperature, where pc ~ {2D)~^ is the density of states for 
conduction electrons with the half bandwidth D. 

Quantum fluctuations should be incorporated for the critical physics at the Kondo 
breakdown QCP, where two kinds of bosonic collective modes will scatter two kinds of 
fermions, that is, conduction electrons and spinous. Gauge fluctuations corresponding 
to phase fluctuations of the hopping parameter Xij = xe*"'-' are introduced to express 
collective spin fluctuations [I5]. Hybridization fluctuations are critical, playing an 
important role for the Kondo breakdown QCP. Such four field variables lead us to 
the following effective field theory in the continuum approximation, 

+ f*{dr - Pf - iar)fa + I (^i " ^^i) f^\'^ 

Znif 

+ b*{dr - /i, - za,)b +-^m- ta,)b\' + ^\b\^ 
+ V{b*clf„ + H.c.) 

+ ^f^^uf^,u + SN{nb + iar), (3) 

where g is an effective coupling constant between matter and gauge fields, and several 
quantities, such as fermion band masses and chemical potentials, are redefined as follows 

A ^ -fit, {2m,)-' = t, {2mf)-' = Jx, 
fic = + '^dt, — fif = ej + \ — 2Jdx- 

Fermion bare bands and e{ for conduction electrons and spinous, respectively, 
are treated in the continuum approximation as €% ~ —2dt + t{k1 + ky + k1) and 
e{ ~ —2Jdx + Jxi^x + + ^^)- '^^^ band dispersion for hybridization can arise from 
high energy fluctuations of conduction electrons and spinous. Actually, the band mass 
of holons is given by m^' ~ NV^Pc/2, where pc is the density of states for conduction 
electrons [HI [12]. Local self-interactions denoted by Ub can be introduced via non- 
universal short-distance-scale physics. Maxwell dynamics for gauge fluctuations appears 
from high energy fluctuations of spinous and holons. 
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Based on the effective Lagrangian, recent studies [HI |T2] developed an Eliashberg 
theory for the Kondo breakdown QCP, where momentum dependence in fermion self- 
energies and vertex corrections are neglected, allowing us to evaluate one loop-level 
quantum corrections fully self-consistently. Actually, this approximation was shown to 
be "exact" in the large limit [I6]. The Eliashberg theory for hybridization fluctuations 
results in the z = ?> Kondo breakdown QCP, discussed in the introduction. 

3. Quantum Boltzman equation study 

We examine electrical transport at the Kondo breakdown QCP based on quantum 
Boltzman equations, where we assume that both hybridization and gauge fluctuations 
are in equilibrium and consider only fermion contributions, consistent with the one-loop 
result for the transport coefficient [1^. Since we have two kinds of fermion excitations, 
we find coupled quantum Boltzman equations for distributions of conduction electrons 
and spinous. Solving such coupled quantum Boltzman equations, we find that the 
diagrammatic result is recovered in the so called "decoupling" limit of these equations, 
where vertex corrections for scattering with hybridization fluctuations can be ignored, 
but those for scattering with gauge fluctuations should be introduced in the spinon 
conductivity. 

Before we perform the quantum Boltzman equation study for the Kondo breakdown 
QCP with two bands, we review on this approach in the U(l) gauge theory with one band 
in order to understand the role of vertex corrections in the transport coefficient p!71 [T8] 
and demonstrate that our treatment successfully recovers the known result [IHl [20] . 

3.1. Application to U(l) gauge theory for a spin liquid state 

We apply the quantum Boltzman equation to the transport problem of U(l) gauge 
theory. 



where D{q,u) = (— ^7^- + Xip'f)" is the gauge propagator with the diamagnetic 
susceptibility Xip and Landau damping coefficient 7^. 

One can obtain this effective field theory from the Hubbard model in the frustrated 
lattice based on the U(l) slave-rotor representation [21], where charge fluctuations are 
gapped at half filling, but magnetic ordering is prohibited owing to the geometrical 
frustration, corresponding to a spin liquid Mott insulator. One also finds this effective 
field theory in the so called algebraic charge liquid for the anomalous normal state of high 
Tc cuprates, derived from the U(l) slave-fermion representation, where spin fluctuations 
described by Schwinger bosons are gapped, but charged excitations represented by 




/ 



dr / d'^r^%ljl{dr - ia^ - + |(<9i - iai)%l)„ 
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fermionic holons are gapless, allowing an anomalous metallic state due to scattering 
with gauge fluctuations [22| [23]. 

Compared with the effective fleld theory for the Kondo breakdown QCP of the 
Anderson lattice model, this U(l) gauge theory is a simplifled version since it does not 
have both holons and conduction electrons. In this section we focus on the mathematical 
structure, in particular, the gauge invariant expression for conductivity [191 120] instead 
of the physical aspect, in order to prepare for the Boltzman equation study of the 
Anderson lattice model. 

We start from the quantum Boltzman equation [23] 

[d^f{uo)]T{k,uj)[A{k,uo)f^^k ■ E = Icoii{k,uj), (5) 

where T{k,uj) and A{k,uj) are the imaginary parts of the retarded self-energy and 
retarded Green's function, respectively, /(co") is the Fermi-Dirac distribution function in 
equilibrium, Vjt is the velocity of fermions, and E is an external electric fleld. Icoii{k, ui) 
is the collision term given by 

Icoii{k,uj) = i:>{k,uj)G<{k,uj) - i:<{k,uj)G>{k,uj), (6) 

where a;) and (^^'^(A;, u;) are lesser and greater self-energies and Green's 

functions, respectively. Using the identity of 

j:>{k,uj)G<{k,uj) - j:<{k,uj)G>{k,uj) 

= 2iT{k, uj)G<{k, u) - iJ:<{k, u)A{k, u), 
where the lesser self-energy is given by 



E<(t,.)=i:r^l^f3Z5(,.. 

„ Jo Tc m,i, 



[{n{u) + l}G<{k + q,u + u)+ n{u)G<{k + q,u- z/)] (7) 

with the Bose-Einstein distribution function n(z/) in the one loop approximation, the 
lesser Green's function is the only unknown function, determined by the quantum 
Boltzman equation. This quantum Boltzman equation is well derived in Ref. |24j . 
based on the Schwinger-Keldysh formulation. 

In the linear response regime we can expand the lesser Green's function up to the 
flrst order of an electric fleld 



G<{k,uj) = iA{k,uj) f{uj) 



E-VkA{k,uj] 



(8) 



^ dco 

where A(fc, u) is the distribution function out of, but near equilibrium due to the electric 
fleld. Inserting this ansatz into the lesser self-energy, we obtain the following expression 
for the lesser self-energy 

du k X q 2 



PC 

^<{k,u) = iY, 



'^D{q,iy)f{uj) 



[[n{u) + f{uj + u)]A{k + q,uj + u)- [n{-u) + f{uj - v)]A{k + q,uj-v)] 
'''^ dv \k X q 2^^^ ^ df{u) 



TT 
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l-/(^ + z/) 



{[n{u) + f{uj + u)] M (A; + q,uj + u)A{k + q,u + u) 

- [n(-z/) + /(a; - u)]^-^^^jf^ A{k + g, ^ - z/)A(fc + g, ^ - z.)}, (9) 



where we used the identities for thermal factors of fermions and bosons, 

{niu) + l}f{uj + u) = f{co){niu) + f{u + z/)}, 
n{u)f{uj -u) = -f{uj){n{-u) + f{uj - u)}, 



and 



Hz/) + i}(- 



df{u + u) 



{n{v) + f{u + v)} 



l-/(^ + z/)/ df{oo) 



1 - /(cj) V duj 
l-f{u;-u). dfioo) 



Inserting both the lesser Green's function and self-energy into the quantum 
Boltzman equation, we find 



A{kF,u) ^ \A{kF,u) + — i 

2 2T{kF,uj) g 



dv 



Hi ^ 



TT 



kF X q 



'^D{q,v) 



[[n{v) + f{u} + v)]A{kF + q,uj + v)- \n{-u) + /(cu - z/)]A(fcj. + g, - z/)} 



A(A;f,u;) 



(10) 



where the momentum is replaced with the Fermi momentum kF because usual transport 
phenomena occur near the Fermi surface except some topological quantities such as Hall 
conductivity |25] and frequency dependence in both the "vertex-distribution" function 
h.{kF-,oj) and thermal Fermi factor is simplified. In this expression the imaginary part 
of the self-energy or scattering rate is defined as 

r-^ dv k x q 2, 



2r{k,u) = Y.f 

„ Jo 



77 



^D{q,iy)\[n{i^) + f{uj + iy)]A{k + q,uj + u) 

(111 



- [n(-z/) + f{uj - u)]A{k + g, - z/)}. 

This approximation will be justified by the fact that it gives rise to the known result in 
the gauge theory context. 

Introducing the relative angle 6 between the initial kF and final kF + q momenta, 
we obtain 

2r{kF,uj) 



A{kF,uj) 



where 



2Ti_cos{kF,uj 
3 



-A{kF,Uj) 



(12) 



2ri_eos(A;F,cu) = j\qq^ f_^d COS e[vlcos\e/ 2)] 



/ —QD{q, u) [1 - cos e] { [n{u) + f{uj + y)]A{kF + g, + z/) 
- [n{-v) + f{uj - v)]A{kF + q,uj ~ v)]. 



(13) 
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In this expression Y^q is replaced with ^fs Jq dqq"^ j\ d cos 6 in = 3, where A is a 
momentum cutoff. 1 — cos^ factor in ri„cos(^F; identifies [2ri_cos(^F; ^)] ""^ with the 
transport time Ttr{uj), capturing large angle scattering dominantly. 

The electrical (charge) or number conductivity is expressed by the lesser Green's 
function, 

Inserting the near-equilibrium ansatz for the lesser Green's function into this expression, 
we obtain the electrical conductivity 

cr,AT) = / ^v,,v,.{-^)A{k,u)A{k,uj) (15) 

where the equilibrium contribution does not generate currents, thus vanishes. 

Inserting the vertex-distribution function into the conductivity expression and 
performing the momentum and energy integration, we reach the final form of the 
conductivity 

a{T) ^ CNFvlrtriT) (16) 

with C = ^ /f^ dy (^2^1)2 ; where Np is the density of states at the Fermi surface and 
the transport time is TtriT) = [2ri_cos(7')]~^, as emphasized before. 

The transport time turns out to be Ttj.{T) oc T~^/^, giving rise to o"(T) oc T~^/^ in 
c? = 3, completely consistent with the previous study [T^l [20] ■ An important point is 
that although the self-energy correction due to gauge fluctuations is diverging at finite 
temperatures, the gauge invariant expression for the conductivity allows only the finite 
result, cancelling the divergence via the vertex correction [T3|. 1 — cos 6* guarantees such 
cancellation. This is the power of the quantum Boltzman equation, imposing the vertex 
correction naturally. 

In this derivation gauge fluctuations are assumed to be in equilibrium. Generally 
speaking, their non-equilibrium distribution due to external fields should be introduced. 
Actually, phonon drag effects are well known in the electron-phonon system |24j . 
Recently, this issue was considered in the spin liquid context with z = 3 gauge 
fluctuations [20], where coupled quantum Boltzmann equations for spinon and photon 
distribution functions are derived. It was argued that such coupled transport equations 
can be decoupled in some cases, where such drag effects are subdominant, compared 
with fermion contributions. 

The present formulation differs from the previous approach in the fact that we 
did not decompose the gauge field as the study of Refs. [TU [20], where the low energy 
gauge field giving rise to divergence is neglected and only high energy gauge fluctuations 
are taken. Although the vertex-distribution function itself is not well defined because 
its part corresponding to the scattering rate is divergent at finite temperatures, we 
found that such decomposition is not necessary because the formal divergence should 
be cancelled in the last gauge invariant physical expression. This spirit goes exactly 
through that of the diagrammatic study. 
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3.2. Application to the Kondo breakdown QCP of the Anderson lattice model 

In the Kondo breakdown scenario we have four kinds of field variables, corresponding to 
conduction electrons, spinous, holons (hybridization fluctuations), and gauge bosons 
(collective spin fluctuations). Our main assumption for the transport study based 
on the quantum Boltzman equation approach is that both hybridization and gauge 
fluctuations are in equilibrium, as pointed out earlier. This assumption is justified 
by the diagrammatic study [1^, where contributions from boson excitations are much 
smaller than fermion contributions, and by the Boltzman equation study of the U(l) 
gauge theory discussed in the previous section. As a result, we are allowed to have two 
coupled quantum Boltzman equations, 

[A,{k, uj)fdJ{uj)Y, ■ v^r,(A;, u) = r^^iiik, cu), 

r^iiik, uj) = 2iT,{k, u;)G< {k, u) - iT.< {k, u;)A,{k, u) (17) 
for conduction electrons and 

[Afik, uj)]'d^fiu;)E . vlTjik, u) = lUk. ^), 

I^k, uj) = 2tTf{k, uj)G<{k, Lo) - zEfik, uj)Aj{k, uj) (18) 

for spinous. 



3.2.1. Contribution of conduction electrons The lesser self-energy for conduction 
electrons arises from scattering with hybridization fluctuations, given by 



j:<{k,u) = v'Y: r-m(g, 

n JO n 



V] 



(19) 



[{n{v) + l}Gj{k + q,uj + v)+ n{u)Gf{k + q,uj- u)] 

in the Eliashberg framework. Since the spinon Green's function appears in the electron 
self-energy, the two quantum Boltzman equations are coupled with each other. This 
coupling effect is the main character for the quantum Boltzman equation of the Anderson 
lattice model at the QCP. 

Inserting the lesser Green's function of spinous 

^dfiuj) 
duj 

into the electron lesser self-energy and the lesser Green's function for conduction 
electrons 

'df{u) 



G<{k, uj) = tAfik, uj) [fiuj) - ■ v{A^(A;, uj) 



(20) 



G<{k,uj) = tA,{k,uj)\f{uj)-(- 



E-vlA,{k,uj) 



^ duj 

into the quantum Boltzman equation for conduction electrons, we obtain 



(21) 



J2 —'^D,{q,u) 



TT 



A.(^^,.).iA.(fc^,.) + ^^,^ ^ 

[[n{u) + f{uj + iy)]Af{k^p, + q,uj + iy) 

- [n{-u) + f{uj - u)]Af{yp + q,uj- z/)}A^(4,u;) 



„,c2 



(22) 
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where 



rOO rill , 

2r,(fc,a;) = V'Y.\ y){\n{v) + /(a; + v)\ 

q Jo TT ^ 

Af{k + q,u; + u) - [n{-u) + f{u; - u)]Af{k + q,uj - u)} 



(23) 



is the scattering rate of conduction electrons and the same approximations as the case of 
the U(l) gauge theory are utilized. It is important to notice that the vertex-distribution 
function for conduction electrons is related with that for spinous. We should know the 
vertex-distribution function for spinous. 



3.2.2. Contribution of spinons The lesser self-energy for spionon excitations results 
from scattering with both hybridization and gauge fluctuations, given by 

r°° dv , 



{k,u;) = E'}<{k,uj) + i:f{k,uj), 



roo fill 

^<{k,u;)^V'Y. -^D,{q,u) 

■' ^ Jo TC 

[{n{u) + l}G<{k + q,u + u) + n{u)Gf{k + q,LJ - p)], 

™^/, X sr^ f°° du k X q 2^ ^ ^ 

A;,u; =E / ' ^Da{q,u) 

■' g Jo TV ruf 



(24) 



[{n{u) + l}Gf{k + q,LO + u)+ n{v)Gj{k + q,uj - u)], 

where the lesser Green's function of conduction electrons appear in the hybridization- 
vertex-induced spinon self-energy while that of spinons arises in the self-energy 
correction via gauge fluctuations. 

Inserting the lesser Green's functions for both conduction electrons and spinons 
into the lesser self-energy and quantum Boltzman equation for spinons, we find 

rf T/2 „oO (111 



2 2r f{k'F,u) q 



{ [n{u) + f{u + u)]Ac{k'p + q,u + u) 

- H-i^) + f{uj - v)]Ac{k% + q,u;- i/)}Ac(A;^, 



+ 



2r/(4, 



dv 



Hi ^ 



TT 



k^p X q 



[[n{v) + f{uj + v)]Af{ki^ + q,uj + v) 

- [n{-v) + /(a; - i^)]A/(A;^ + a; - y)}^^f{kU ^), 



(25) 



where 



2Tf{k, uj) = V^Y. —^D,{q, u){[n{u) + f{uj + u)]A,{k + g, + i/) 

n Jo n 



+ E 



H-T^) + f{u;- u)]Ac{k + g, a; - I/) } 
^ dz/ k X q 2 



TT 



^Da{q, i^){[n{u) + f{uj + p)]Af{k + q,uj + v) 
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- [n{-v) + f{u-v)\Af{k + q,uj-v)} 

= 2T){k,uj) + 2T){k,uj) (26) 

is the scattering rate of spinons resulting from scattering with both hybridization 
T^f{k,ijj) and gauge fluctuations r^(A;,a;). One can check Eqs. (25) and (26), considering 
that the hybridization-induced part is basically the same as that of the quantum 
Boltzman equation for conduction electrons and the gauge-fluctuation part coincides 
with that shown in the U(l) gauge theory of the previous section. 

Inserting the vertex-distribution function [Eq. (22)] for conduction electrons into 
the vert ex- distribution function [Eq. (25)] for spinons, we find the following expression 
for spinons 

KfikUu^) = \{Af{kUu^) + ^-J^^f^^A^{k%,u^)} 



with 



+ 1 J. lU ^ V ih'^ ,{\ ^ r (uf V)^^^' > ^ > 

roo , 

/ —QDb{q, [n{u) + f{uj + v)]A^{k% + q,uj + v) 

Jo TT 

-[n{-iy) + fiu- iy)]A^{k''p + q,uj-iy)], 

2r/,cos(4, ^) = 2A3 io ^^^^ J-i ^ ^^-^ 

/ cos2(^^^^/2)]5>D„(g, v){[n{v) + f{uj + z/)] 

Af{yp + q,uj + v) - [n{-v) ^ f{u - v)]Af{yp ^ q,u - v)], 

2rc,cos(^F>'^) = dqq^ J ^d COS Ocfi^^ COS Ocf) 

r°° du ( f 

/ — m(?, [n{iy) + f{uj + v)\Af{k^p + + I/) 

Jo 71 

- [n{-u) + f{uj - u)]Af{ki + q,uj-u)}, (28) 

where 6cf represents an angle between the initial electron velocity Vp and final spinon 
velocity Vp and Off is defined in the similar way, but between spinons. 
We obtain the spinon vertex-distribution function 

A/(4,c^) = l{rf{k{.,u)Af{kfp,u) + r'f^,Jkfp,u;)A,{k'p,u^ 

{r}(4,c^) + ri,_^Jk{.,u;) - rJ,eos(4,^)%7^%^}"\ (29) 



where 



3 , o 



2r^,i-cos(^F, ^) = ^j^ '^11^ J_^dcos9ff[l - COS 9 ff] 
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TT 



/p^ cos\d ff/2)]^Da{q, iy){[n{u) + f{u + v)]Aj{k^p + q,u + v) 
- [n{-u) + f{uj - u)]Af{yp + g, ^ _ z.)} (30) 
is identified with [rj'J.(ct;)]^^ as shown in the U(l) gauge theory of the previous section. 

3.2.3. Conductivity in the decoupling limit In the vertex-distribution function for 
spinous [Eq. (25)] we neglect the couphng term A.^k'^p^uj) as the zeroth order 
approximation for the transport study, named as the decouphng hmit. One may 
understand vahdity of this approximation, based on the fact that spinous are heavily 
massive denoted by a ^ 1 and scattering with conduction electrons will not affect their 
dynamics much. Then, we find 



1 Vf{kfp,u)A^{y 



Inserting this expression into the spinon conductivity, we obtain 

CN^pv{,^ 

^^^^^ - 2r^(T) + 2r«,,_(T)' ^^^^ 

exactly the same as that of the diagrammatic study [lAl, where vertex corrections are 
introduced only for the scattering channel with gauge fluctuations. As discussed in 
the previous section, z = 3 gauge fluctuations give rise to divergence for self-energy 
corrections to spinous, but cancelled by vertex corrections, allowing the gauge invariant 
finite physical conductivity proportional to ~ T"^/^ in d = 3. One may ask why the 
same situation does not happen for scattering with z = 3 hybridization fluctuations. 
Actually, such 2; = 3 dynamics of holons is cut by an intrinsic energy scale £"*, and 
scattering with z = 2 hybridization fluctuations below E* does not cause the divergence 
for self-energy corrections. 

The vertex-distribution function for conduction electrons becomes 

2 2 T,{k%u;) V){yp,u) + Tl^_,,M,u) 

where scattering with spinous is incorporated through the vertex-distribution function 
for spinous because hght conduction electrons can be much affected. However, calling 

= - « « 1, (34) 

the second contribution in the electron vertex-distribution can be neglected. As a result, 
the conductivity from conduction electrons is free from vertex corrections, becoming 

a^iT) = C^f^, (35) 

^ ^ 2r,(T)' ^ ^ 

which coincides with that of the diagrammatic study ^l4j showing Tc{T) ~ Tln(T/£'*) 
in the z = 3 critical regime. 
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The last work is to find an actual expression for the physical conductivity, referred 
as the loffe-Larkin composition rule [26] 

where (Tb{T) is the ho Ion conductivity, much smaller than fermion contributions 
justifying the last approximation. 

One can ask the role of the spinon conductivity for any physical response functions. 
Actually, it contributes to the physical thermal conductivity given by the corresponding 
loffe-Larkin composition rule 

^^^^K^_ (37) 

where Kcj(T) are thermal conductivity of conduction electrons and spinous, respectively, 
and holon contributions are also neglected. Assuming that the Wiedemann-Franz law 
holds for each fermion sector, proven to be correct at least in the one loop approximation 
[Ti] , we find 

^-yMT)+a,(T)), (38) 

suggesting that the Wiedemann-Franz law should be violated due to the presence of 
additional entropy carriers, that is, spinous at the Kondo breakdown QCP in the low 
temperature limit, i.e.. 



with Lq = the value of the Fermi liquid. This result would be robust beyond our 

approximation because this expression includes just density of states and velocity at the 
Fermi energy, thus expected to be governed by a conservation law. 



4. Summary 

In this paper we developed the quantum Boltzman equation approach for the Kondo 
breakdown QCP, which two bands for conduction electrons and localized fermions are 
involved with, where scattering with z = 3 critical hybridization fluctuations and ^ = 3 
gapless gauge bosons relaxes their dynamics. Our main problem was to understand 
the role of vertex corrections in their transport phenomena, crucial for the T-linear 
non-Fermi liquid resistivity in the z = 3 critical theory. 

Only one band of spinous is involved for scattering with gauge fluctuations, and 
their associated vertex corrections are introduced in the usual way, demonstrated in 
the U(l) gauge theory of section 3-1. Our treatment for gauge fluctuations is different 
from the previous study [19] , in respect that the vertex-distribution function is not well 
defined owing to its formal divergence associated with vertex corrections, but cancelled 
in the physical conductivity through self-energy corrections as it should be, consistent 
with the diagrammatic approach [18], while divergent contributions are thrown away 
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in the vertex-distribution function of the previous study, well defined. Of course, both 
approaches give the same result. 

On the other hand, two bands should be taken into account for scattering with 
hybridization excitations, giving rise to coupled quantum Boltzman equations. In the so 
called decoupling limit where coupling effects are neglected for the Boltzman equation of 
spinons while they are allowed for that of conduction electrons, the vertex correction for 
conduction electrons associated with hybridization fluctuations turns out to be irrelevant 
due to heavy mass of spinons denoted by o; ^ 1. Results of the diagrammatic approach 
are recovered from the quantum Boltzman equation approach in the decoupling limit. 

The next task is what happens beyond the decoupling limit. Our preliminary 
analysis shows that vertex corrections seem to appear in the scattering channel with 
hybridization fluctuations. However, we do not find the corresponding diagram for such 
a correction at present. It remains as an interesting future study. 
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